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It is known that the coordination number (CN) of atoms or ions in many materials increases through 
application of sufficiently high pressure. This also applies to glassy materials. In boron- containing glasses, 
trigonal B0 3 units can be transformed into tetrahedral B0 4 under pressure. However, one of the key 
questions is whether the pressure-quenched CN change in glass is reversible upon annealing below the 
ambient glass transition temperature (T g ). Here we address this issue by performing n B NMR 
measurements on a soda lime borate glass that has been pressure-quenched at —0.6 GPa near T g . The results 
show a remarkable phenomenon, i.e., upon annealing at 0.9 T g the pressure-induced change in CN remains 
unchanged, while the pressurised values of macroscopic properties such as density, refractive index, and 
hardness are relaxing. This suggests that the pressure-induced changes in macroscopic properties of soda 
lime borate glasses compressed up to —0.6 GPa are not attributed to changes in the short-range order in the 
glass, but rather to changes in overall atomic packing density and medium-range structures. 

Any natural or man-made material will exhibit a significant change in microstructure and properties when 
subjected to high pressure 1 " 5 ; particularly the coordination number (CN) of atoms or ions will generally 
increase with pressure 6 " 9 . It has long been known that when a glassy material or its liquid state is subjected 
to sufficiently high pressure, significant changes can take place in the local and medium- range structure, vibra- 
tional density of states, and physical properties 4 " 81011 . In recent years, theoretical modeling approaches have 
provided insights regarding the structural response of various non- crystalline materials to composition, temper- 
ature, and pressure 11 " 15 . Nevertheless, the relationship between structure and dynamics in glass and glass-forming 
melts under high pressures remains a challenging problem in condensed matter science. The difficulties in 
addressing this problem are due to both experimental limitations at high pressures and the inadequacy of 
computer simulations under such conditions. It is crucial to determine the link between the microscopic structure 
and macroscopic properties of glasses under high pressure from both scientific and technological perspectives, 
since the glass structures frozen-in under elevated pressure may give rise to properties unattainable under 
ambient pressure. Borate glass is the quintessential example of glass capable of dramatic changes in short-range 
order as a function of composition, pressure, and thermal history, and hence we select a simple ternary borate 
system, namely, sodium- calcium-borate glass, for this study in order to address the above mentioned challenging 
problem. Boric oxide (B 2 0 3 ) is also widely used as a network forming constituent in many high-tech glass 
materials owing to its contribution to high glass forming ability and low melting temperature, and for its favorable 
impact on thermal, mechanical, and optical properties. 

One of the most striking structural features of borate glasses is the transformation of the coordination number 
(CN) of boron from three to four upon compression 16 " 19 , which is associated with anomalous pressure depend- 
ence of viscosity 616 and topological disorder 812 . Densification of borate glasses can be achieved by applying 
isostatic pressure to glass directly at room temperature 1617 , or by pressure-quenching from the molten state 18 " 21 . 
However, the elastic part of the densification relaxes during pressure release, i.e., the density of a compressed glass 
at ambient pressure decreases to the value before compression upon reheating below the glass transition tem- 
perature (T g ) 17 . For vitreous B 2 0 3 , it has been shown that the fraction of boroxol rings decreases with increasing 
pressure 16 ' 20 ' 22 , leading to densification of the glass 16 and increase of elastic moduli 20 . When the pressure applied at 
room temperature is released, such densification is irreversible since the local structures may be topologically/ 
stereochemically unfavorable for the reformation of boroxol rings 16 . However, it is not clear whether this irre- 
versibility of densification is linked to the CN change of boron 23 . More recently, it has been found that the B0 3 
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fraction suddenly drops at pressures just above 4 GPa and then 
approaches zero as the pressure is further increased 17 . A structural 
response of vitreous B 2 0 3 to pressure below 4 GPa was not detected. 
After decompression from high pressure the boron coordination 
reverts back from tetrahedral to trigonal, however, the densification 
is apparently permanent 17 . It has been attempted to correlate these 
changes in microscopic structure with the macroscopic properties. 
For example, the viscosity of the B 2 0 3 liquid along the melting curve 
has been shown to decrease by 4 orders of magnitude upon a pressure 
increase up to 5.5 GPa but then remains unchanged upon further 
increase of the pressure 22 . However, a generally accepted viewpoint 
about the microscopic origin of the pressure-induced changes in the 
macroscopic properties is still lacking. 

The abovementioned studies were carried out on vitreous B 2 0 3 at 
room temperature and rather high pressure up to 22.5 GPa. It should 
be noted that the structure of glasses compressed at room temper- 
ature (in diamond anvil cell) is different from that obtained by pres- 
sure-quenching liquids from above T g 21 , even though the main 
structural changes upon compression (i.e., decreased fraction of bor- 
oxol rings and increased fraction of tetrahedral boron) are identical. 
It should also be noted that modified borate and pure B 2 0 3 glasses 
have different initial concentration of boroxol rings and may exhibit 
different densification mechanisms. However, Lee et al. have shown 
that the pressure dependent structural changes of a lithium borate 
glass are similar to those of vitreous B 2 0 3 24 . In particular, they 
showed that the pressure-induced CN change in Li 2 B 4 0 7 glass from 
three- coordinated to four- coordinated boron at room temperature 
begins around 5 GPa and the B0 4 fraction increases with pressure 
from about 50% at 1 atm to more than 95% at 30 GPa. 

In this work we investigate the microscopic and macroscopic res- 
ponses of the soda-lime borate glass to comparatively low pressure at 
an elevated temperature around the T g , since numerous high pres- 
sure studies at room temperature have already been reported in 
literature and it is also our aim to detect the sensitivity of glass 
structure and properties to the low pressure. Moreover, these condi- 
tions are chosen to be able to prepare compressed samples, which are 
large enough to allow for subsequent characterization of mac- 
roscopic properties. The as-produced glass has a nonzero equilib- 
rium concentration of B0 4 due to the presence of network modifiers 
(Na 2 0 and CaO) that act to charge balance the tetrahedral boron 
species. Only recently has it been found that boron speciation (i.e. the 
CN) changes occur at rather low pressure (<0.6 GPa) near T g and 
can be frozen-in under pressure 18 . Here the key questions arise: Is the 
pressure-quenched CN change reversible upon annealing at tem- 
peratures slightly below T g at atmospheric pressure? Are the pres- 
sure-induced changes in the macroscopic properties also reversible 
upon annealing? What is the microscopic origin of the pressure- 
induced changes in the macroscopic properties? We answer these 
questions by investigating the structure and property responses of 
the sodium-calcium-borate glass to pressurization at <0.6 GPa and 
subsequent annealing at 0.9 T g under ambient pressure. We thereby 
hope to obtain a better understanding of the pressure-structure- 
property relationship of this and other glass systems, which can be 
used for tailoring both microscopic and macroscopic properties of 
glassy materials. 

Results 

Physical properties and heat capacity. The glass composition under 
study is 25Na 2 0 - lOCaO - 65B 2 0 3 (mol%). As reported 
previously 25 , the density and Vickers hardness of this glass increase 
approximately linearly with increasing isostatic pressure (inset of 
Fig. 1). The glass isostatically compressed at 570 MPa is then 
annealed under ambient pressure at 688 K, i.e., 0.9 times its 
ambient T g for various durations (f a ). This results in a decrease of 
both density and Vickers hardness with increasing annealing 
duration (Fig. 1). In other words, these macroscopic properties of 
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Figure 1 | Relaxation of macroscopic properties. Evolution of density (p) 
and Vickers hardness (H v ) with annealing duration (f a ) at 0.9 T g = 688 K 
of the borate glass compressed at 570 MPa. The dashed lines indicate the 
values of p and H v prior to compression. Inset: impact of isostatic pressure 
(p) on p and H v . The errors of p and H v are around ±0.005 g/cm 3 and 
±0.2 GPa, respectively. 

the compressed glass are relaxing during annealing towards those of 
the glass prior to compression. While Vickers hardness has 
essentially decreased to its original value prior to compression 
(5.0 GPa) after annealing for 6 h, density has only decreased to 
—2.47 g/cm 3 after annealing for 24 h, which is significantly higher 
than the value prior to compression (2.438 g/cm 3 ), i.e., density has 
relaxed by only 43% of the total possible relaxation. 

With increasing isostatic pressure, the overshoot in the isobaric 
heat capacity (C p ) above T g increases as illustrated in Fig. 2a, where 
the evolution of C p with temperature during the first DSC upscan is 
shown. This overshoot is considered to be a direct consequence of the 
nonexponentiality of the relaxation process 26 , i.e., due to broadening 
of the relaxation time distribution. This also implies that subjecting 
the glass to high pressure forces the glass into a lower region of the 
enthalpy landscape compared to the glass under ambient pressure, 
even though it should be noted that the enthalpy landscape itself is 
changing as a function of pressure. In agreement with previous stud- 
ies 19,27 , we find that an increase in isostatic pressure enhances the 
nonexponentiality of the enthalpy relaxation as measured ex situ. 
With increasing pressure on the glass, the density of the glass 
increases, and at the same time the topological degree of atomic 
freedom decreases as well due to the B0 3 ^B0 4 conversion, i.e., 
the increase of network connectivity. Upon heating through the glass 
transition region, the glass compressed at higher pressures should 
exhibit a larger jump in configurational entropy in order to approach 
the liquid state. This is also shown in Fig. 2b, in which the pressure 
dependence of the fictive temperature (Tf) is plotted. The procedure 
for determining T f based on DSC data is described elsewhere 28,29 . 

The compressed glass relaxes in the glass transition region during 
the first upscan and recovers back to its original state with respect to 
enthalpy, i.e., that of the uncompressed glass. In other words, the 
enthalpy level of the glass reaches that of the glass cooled under 
standard conditions, e.g., at 10 K/min and ambient pressure. This 
is seen from the second DSC upscans as shown in the inset of Fig. 2a 
and in Supplementary Figure SI. Enthalpy recovery of the com- 
pressed glass during the DSC upscan is a "structural depression" 
process, induced by thermal excitation under ambient pressure 27 . 
The structural depression leads to an increase in enthalpy, which is 
manifested as an enhancement of the overshoot during the first DSC 
upscan. This is also evident from the relaxation of refractive index 
before and after the first DSC upscan (Fig. 3). With increasing iso- 
static pressure, the refractive index at 633 nm increases due to den- 
sification, but the measured value on the sample following the first 
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Figure 2 | Pressure history dependence of calorimetric glass transition. 

(a) Evolution of endothermic overshoot in isobaric heat capacity (C p ) at 
the glass transition for samples with different pressure history. DSC scans 
were performed at 10 K/min at ambient pressure. Inset: second DSC 
upscans at 10 K/min of all samples following prior up- and downscans at 
10 K/min. (b) Pressure history dependence of the Active temperature (T f ) 
evaluated from the first and second DSC upscan, respectively. 

DSC upscan is independent of the initially applied pressure during 
compression. 

Structural response. n B magic-angle spinning (MAS) NMR spectra 
obtained at 16.4 T (700 MHz) for the 0.57 GPa compressed glass 
annealed for different durations are shown in Fig. 4a. These spectra 
are characterised by a broad peak centered at +15 ppm, 
corresponding to B m sites, and a relatively narrow peak centered 
around +2 ppm, corresponding to B IV sites. The spectra vary 
slightly with annealing time, which reflects minor changes in either 
the relative proportions of B m and B IV , or changes in bond angles and 
distances involving boron and oxygen. We quantify these differences 
by accurate simulation of the spectral lineshapes and subsequent 
determination of the fraction of tetrahedral to total boron (N 4 ) 
through integration. The simulation parameters are given in 
Supplementary Table SI, and two examples of the deconvolution 
are shown in Supplementary Figure S2. Prior to annealing, we find 
that with increasing isostatic pressure, N 4 increases from 44.2 to 46.5 
at% (inset of Fig. 4b and Supplementary Figure S3a). N 4 also 
increases upon compression in pure B 2 0 3 glass 17 and in boron- 
containing multicomponent E-glass 18 . In Fig. 4b we observe an 
interesting phenomenon, viz., N 4 of the compressed sample 
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Figure 3 | Recovery of refractive index during DSC scan. Pressure 
history dependence of the refractive index (n) at 633 nm before and after 
the first DSC scan, as shown in Fig. 2a. 

remains almost constant within the error range when extending 
the annealing time up to 25 hours at the annealing temperature of 
688 K (0.9 T g ) and at the ambient pressure. In other words, the glass 
structure, as defined by boron coordination number, does not return 
to its original state prior to compression upon annealing. Although 
N 4 does not appear to change, the u B MAS NMR spectra do show 
differences with annealing time. For example, the base of the 
tetrahedral boron resonances in Fig. 4a exhibit some annealing 
time dependence, signifying a change in environment of the B IV 
groups, possibly due to changes in bond angle distribution around 
these boron atoms. The B m features in Fig. 4a do not show much 
change with annealing time, indicating that annealing has little 
impact on B-O bond distance or angles for these structural elements. 

Supplementary Figure S3b shows 23 Na MAS NMR spectra of the 
samples prior to annealing. We see a small but systematic increase in 
frequency shift with increasing pressure, which is due to a decrease in 
the mean Na-O bond distance upon compression 18 . Fig. 4c shows the 
23 Na MAS NMR spectra for the compressed- annealed glasses. These 
MAS NMR spectra overlap completely, and to further study the 
sodium speciation, 23 Na triple-quantum magic-angle spinning 
(3QMAS) NMR data were also collected. Isotropic projections from 
these data (Supplementary Figure S4a) show very little change with 
annealing time, and thus, the estimated quadrupolar coupling prod- 
uct (P Q ) and isotropic chemical shifts (5 iso ) for 23 Na also do not 
indicate any detectable change in sodium environment in the com- 
pressed glass upon annealing (Supplementary Figure S4b). This sug- 
gests that the Na-O bond length does not change with annealing 
time, and hence it does not determine the recovery of the mac- 
roscopic properties upon annealing. 

Discussion 

Our work has shown that upon annealing at 0.9 T g of a soda lime 
borate glass compressed up to —0.6 GPa, the pressure-induced 
change in CN remains unchanged while the pressurised values of 
macroscopic properties such as density, refractive index, and hard- 
ness are relaxing. This suggests that the pressure-induced changes in 
macroscopic properties of such glasses are not attributed to changes 
in the short-range order, in agreement with previous studies at rela- 
tively low pressure 30 . In sodium boroaluminosilicate glasses it has 
been shown that while pressure-induced density changes are related 
to changes in boron coordination, the relatively small difference in 
partial molar volume of the B0 4 and B0 3 structural units cannot 
account for the dominant part of the density change 18 . Instead it was 
suggested that the shortening of the Na-O bond upon compression is 
a more likely factor governing the density change 18 ' 31 . Such changes 
were observed at considerably higher pressures than the current 
study and thus may not in fact account for the findings reported 
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Figure 4 | Structural response to annealing, (a) Solid state n B MAS NMR 
spectra at 16.4 T of the samples compressed at 570 MPa (solid line) and 
then annealed for various durations (t a ) at 0.9 T g = 688 K. (b) Evolution of 
fraction of tetrahedral to total boron (N 4 ), which is obtained through 
integration of peaks in Fig. 4a, with annealing duration. Inset: impact of 
isostatic pressure on N 4 in absence of annealing (corresponding U B spectra 
are shown in Supplementary Figure 2a) . The error of N 4 is ± 0.5 at% for the 
variable pressure series (inset) and ± 1 at% for the variable annealing time 
study. The latter has higher error due to changes in lineshape and thus 
more difficulty in fitting the n B MAS NMR data, (c) Solid state 23 Na MAS 
NMR spectra at 16.4 T of the samples compressed at 570 MPa and then 
annealed for various durations at 0.9 T„ = 688 K. 



herein. For example, Allwardt et al reported a change in 23 Na chem- 
ical shift of 3.5 ppm at 10 GPa, roughly corresponding to a 0.005 nm 
decrease in Na-O distance 32 . The chemical shifts for sodium in these 
soda-lime borates vary slightly with annealing, but within roughly 1- 
2 ppm and due to measurement uncertainty, appear to indicate a 
negligible contribution of Na-O bond distance changes to the density 
change. Shannon & Prewitt 33 indicate that Na-O bonds should be 
more compressible than Ca-O bonds. Since Na-O bonds do not 
change in the detectable range, the contribution of Ca-O bond con- 
traction under pressure could thus likely be ruled out, and therefore 
the expansion of these bonds does not occur upon annealing, i.e., 
should not be the reason for the density decrease. 

Although boron speciation and modifier-oxygen bond lengths do 
not change with increased annealing time, density partially 
decreases. This might be explained as follows. In addition to coordi- 
nation numbers and bond distances, bond angle distributions (B-O- 
B and O-B-O) and superstructures (e.g., boroxol rings 34 ) may change 
upon annealing 16 . The change in bond angle and configuration of B IV 
groups with annealing time could be related to the evolution of the 
B IV resonance around the frequency shift of +2 ppm in the n B NMR 
spectra (Fig. 4a). By isostatic compression, not only the connections 
of neighbored B0 4 tetrahedra, but possibly also the linkages between 
boroxol rings, become stronger or tighter, and this is accompanied by 
the narrowing down of the bond angle distribution. Consequently, 
the total potential energy in the glass decreases with increasing the 
compression, which is why the glass transition overshoot becomes 
larger upon DSC upscanning (see Fig. 2a). Once the compressed glass 
is subjected to annealing around T g , the "tightened state" of glass will 
relax, leading to the decrease of density and hardness. However, the 
thermal energy at 0.9 T g could be insufficient for breaking the B-O 
bonds in favor of formation of B0 3 units, despite the fact that struc- 
tural relaxation generally occurs even at temperatures below 0.9 T g 35 . 
It should be noted that changes in boron speciation due to annealing 
at 0.92 T g have previously been reported for non-compressed boro- 
silicate glasses 36 and the thermomechanically trapped state of the 
glass studied herein thus appears to be different from that of the 
thermally trapped glass. 

It is known that the a relaxation is decoupled from the ft and y 
relaxation below T g 37 . Generally, for network glass systems, the a 
relaxation is dominated by the change of the network connectivity 
(e.g., the boron speciation), and the ft relaxation is controlled by the 
local motion of structural units (e.g., the change of bond angles). In 
this context, it is understandable that the density relaxes upon 
annealing slightly below T g , whereas the CN remains unchanged 
since the former is of the relaxation feature requiring low temper- 
ature and the latter is of the a relaxation feature demanding the 
temperature to exceed T g . According to literature the B m -to-B IV ratio 
increases with increasing temperature above T g 36 . Thus, it is expected 
that the thermomechanically induced B IV units would be converted 
to B 111 units when the temperature is well above T g . In contrast to the 
density decay, the hardness is fully recovered by annealing at 0.9 T g to 
the original value more quickly, as shown in Fig. 1. This could be 
related to the dominant contribution of the y relaxation process, or to 
contributions of hydrated structure relaxations of a near surface 
volume (OH-groups and molecular water) to the relaxation of hard- 
ness. The latter has been found recently to proceed more rapidly than 
both a and /? relaxations 38 . This implies that the hardness relaxation 
be governed by the fast local network relaxation of the surface layer 39 . 

The irreversibility of the thermomechanically induced change in 
boron speciation upon annealing at 0.9 T g implies that the free energy 
barrier for forming a given structural state (i.e., boron speciation) can 
be overcome not only chemically, but also thermomechanically. This 
allows for tailoring of both microscopic and macroscopic properties, 
since the increased network connectivity that remains after anneal- 
ing should increase, e.g., chemical durability and thermal shock 
resistance. To further understand this effect, we consider the 
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enthalpy landscape view of glasses 40 . At high temperatures (i.e., well 
above T g ), the system can flow freely among its configurational 
microstates, corresponding to the case of an ergodic, equilibrium 
liquid. As the system is cooled, the configurational transitions occur 
less frequently owing to the loss of thermal energy. At the glass 
transition, there is a continuous breakdown of ergodicity as the sys- 
tem gradually becomes trapped in a subset of the available config- 
urational phase space known as a "metabasin", i.e., a group of 
configurations that are mutually accessible at a given temperature 
and for a given observation time 41 . The relatively slow configura- 
tional transitions among different metabasins typically involve a 
concerted series of such transitions. During annealing of the borate 
glass in this work, there is more thermal energy available for the 
system to sample the phase space. However, since there is no decrease 
of N 4 during annealing, the boron coordination state is thermome- 
chanically trapped in a basin in the energy landscape as a con- 
sequence of the isostatic compression although the macroscopic 
properties relax on the time scale of the annealing time. This kinetic 
effect is caused by a high activation barrier in the enthalpy landscape, 
preventing the boron coordination number from decreasing. 
Another possibility is the lack of a thermodynamic driving force 
for the coordination change. According to the random pair model 
of Gupta, the equilibrium N 4 value of the studied composition is 49.2 
at% 42 ' 43 . Hence, the N 4 of the compressed glass is closer to the ambi- 
ent equilibrium value than that of the as-prepared glass (see inset of 
Fig. 4b), i.e., there is no thermodynamic driving force to convert B IV 
into B m units. 

Methods 

Sample preparation. Glass with composition (in mol%) of 25Na 2 0 - lOCaO - 
65B 2 0 3 was prepared using melt-quenching technique in an inductively heated 
furnace, as described in details elsewhere 25,43 . Sample rods (4 X 4 X 30 mm 3 ) were 
then isostatically compressed using a cold seal pressure vessel under argon gas. For 
details of the setup, see Ref. 25. The samples were heated under pressure to a 
temperature around T g + 20 K (T g = 764 K), followed by equilibration at this 
temperature for ~3 min, and finally cooled to room temperature at an initial cooling 
rate of ~3 K/min. The experiments were performed at the following pressures (p): 
0.1, 100, 200, 300, 400, 500, and 570 MPa. Relaxation studies were performed on the 
sample compressed at 570 MPa. This was done by isothermal heat-treatment at 0.9 T g 
= 688 K for the following durations (f a ): 15, 30, 60, 120, 180, 360, and 1440 min. 

Physical property measurements. Density (p) was measured in ethanol using 
Archimedes' principle. Vickers hardness (H v ) was determined using a micro- 
indentor (HMV2000, Shimadzu) operated at a load of 9.81 N at ambient conditions 
using a dwell time of 15 s. The diagonals of the indentation were measured using a 3D 
laser scanning microscope (VK-9700K, Keyence). Refractive index (n) of samples 
before and after the DSC measurements (see below) were performed at 633 nm using 
a low range Precision Refracto meter (Bausch & Lomb). 

Differential scanning calorimetry. Isobaric heat capacities (C p ) of the investigated 
samples were determined using a differential scanning calorimeter (DSC Netzsch 
404C). DSC runs included measurements of the baseline (two empty Pt-Rh crucibles, 
6 mm in diameter, 0.1 mm wall thickness covered with a lid), a sapphire standard 
(with one crucible containing the standard and the other empty), and finally the 
sample (with one crucible containing the sample and the other empty). The glass 
samples were polished to within 1 urn to ensure an accurate fit with the bottom of the 
crucible and to reach a mass comparable to that of the sapphire standard (—56 mg). 
Calorimetry was performed under argon flow of 40 ml/min and the glasses were 
placed on the DSC sample holder at room temperature. Subsequently they were 
heated to an initial temperature of 313 K, held isothermally for 15 min, and then 
heated further at a rate of 10 K/min to a temperature 60 K above the T g in order to 
fully relax the sample. At the next step, the sample was cooled at a rate of 10 K/min to 
313 K, and then held for one hour prior to the following heating, which occurred at a 
rate of 10 K/min. 

Solid state NMR. U B and 23 Na MAS NMR experiments were conducted at 16.T using 
a commercial spectrometer and MAS NMR probes. Resonance frequencies for U B 
and 23 Na at this external magnetic field strength were 224.51 and 185.10 MHz, 
respectively. Samples were crushed using an agate mortar and pestle, packed into 
3.2 mm zirconia rotors and spun at frequencies of nominally 20 kHZ. U B and 23 Na 
MAS NMR spectra were collected using short radio -frequency pulses (0.6 |is, 
equivalent to tt/12 tip angles), relaxation delays of 2 seconds and signal averaging of 
1000 to 2000 acquisitions. Data were frequency referenced to aqueous boric acid at 
19.6 ppm and aqueous NaCl at 0 ppm for n B and 23 Na, respectively. n B MAS NMR 
spectra were fit using DMfit 44 to reproduce B m and B IV lineshapes, and N 4 values were 



determined from integration of these resonances and with consideration of a small 
correction in B IV intensity due to overlapping satellite transitions 45 . 23 Na 3QMAS 
experiments were conducted also at 16.4 T (185.10 MHz resonance frequency) using 
a 3.2 mm MAS NMR probe and sample spinning of 20 kHz. The pulse sequence used 
two hard rf pulses (2.3 and 0.9 us), followed by a z-filter echo using a delay of 10 |is 
and a soft reading pulse of 15 us. 
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